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Abstract 


Slicing a signal’s spectrum into channels is a common requirement in telecommunications and signal processing. Polyphase 
FFT Channelizers are digital structures that are very efficient at this task. Consequently, they are widely used in various 
areas of digital signal processing. Since their development in 1974, detailed derivations have been scarce and existing ones 
are often very compact or focus on optimized versions of the channelizer, which are harder to grasp. This paper details a 
new derivation for the non-optimized, basic channelizer as well as the optimized version with a commutator and one with 
a double commutator for reduced aliasing. The reader requires only basic system theory knowledge. 


1 Introduction 


Polyphase FFT Channelizers, first described in 1974 
by Bellanger and Daguet [1], require far less calcula- 
tion effort channelizing signals than a cascade of sin- 
gle filters could [2]. They can construct channelizers 
for television [3], receivers for FM [4] and OFDM sig- 
nals [5] and have been suggested as key components 
for generic wideband receivers [6, 7]. Furthermore, 
they are also a valuable asset for wavelet analysis [2]. 

Unfortunately, complete and detailed derivations 
for these complicated structures have been scarce. 
For instance, research papers [8] and [6] discuss op- 
timized versions of the channelizer, with the rate re- 
duction being the first block in the flowgraph. Nat- 
urally, such optimized versions are harder to under- 
stand. Conversely, this paper presents both funda- 
mentals and optimizations. Section 2 introduces the 
necessary mathematical foundations and system the- 
ory knowledge. Afterwards, Section 3 explains chan- 
nelizers’ fundamentals and describes how a reused 
(modulated) model filter can construct a filter bank; 
Section 4 provides a detailed explanation of how 
polyphase decomposition of said model filter con- 
structs the Polyphase FFT Channelizer. Finally, Sec- 
tion 5 shows how the resulting non-optimized struc- 
ture can be optimized by moving the sample rate 
reduction step to the start of the system’s process- 
ing chain. Additionally, Section 5 establishes foun- 
dational knowledge necessary to implement further 
optimizations and introduces one such highly opti- 
mized version, a channelizer with a reduced sam- 
ple rate reduction and, consequently, reduced self- 
aliasing. 
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2 Mathematical Background 


This section shall briefly establish the necessary 
mathematical and system theory knowledge back- 
ground. 


2.1 Conventions 


All signals in this work are represented as com- 
plex numbers, while all filter taps are real-valued 
if not stated otherwise. Functions, variables and se- 
quences are not marked explicitly to be complex. 

Throughout this paper the following mathemati- 
cal markers are used: 


¢ Exponent (p), ie., FY), identifies a function to be 
a polyphase component of the original function 
F, 

¢ The caret symbol, ie., X,, denotes signals as fil- 
tered by a polyphase-filter, without themselves 
being polyphases. 

¢ The tilde, ie., 7, marks temporary functions or 
variables that are created for formal reasons. 


Frequencies, including the sampling frequency Og 
and the Nyquist Frequency Ow, are presented as nor- 
malized frequencies  [9, 2, 10], allowing for the for- 
mulation of equations independent of sampling time 
Ts: 


O = wTs (1) 
with Os = wsTs = 27 and ON = 1 (2) 


When sample rate changes are discussed, how- 
ever, fs is preferred because a normalized frequency 
would be less intuitive in this context. 
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(a) Filtering with subsequent reduction. 
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(b) Reduction with subsequent filtering. 


Figure 1: Equivalence of two positions of rate reductions. The 
signals X and Y are identical for both versions (a) and (b), but 
the filter changes. 


2.2 Sample Rate Changes 


It is preferable to keep sample rates fs as low as pos- 
sible in digital systems because electronics used in 
modern digital circuits consume more energy with 
increasing clock frequency [9, 2]. Reducing fs be- 
fore a signal processing block such as a filter is espe- 
cially useful as it reduces the number of operations 
per second, resulting in a far more efficient system 
structure. This is possible thanks to the Noble Iden- 
tity, which defines several equivalences of sample 
rate change positions [2, 11], some of which are pre- 
sented here. 

A conventional approach to rate reduction is re- 
ducing the rate at the exit of a system, as shown in 
Fig. la. However, it is possible to reduce fg before 
a signal processing block by using the equivalence 
from Fig. 1b. As shown, moving the rate reduction 
requires that the system’s transfer function is either 
expanded or reduced by M, depending on which 
system is the reference system. The relationship be- 
tween sample rate changes in time and frequency 
domain is given by [2, 11, 10]: 


F(z) eo f{[n| ©) 
F (zi) eo f[m] = f[nM] (4) 
AER ee 


Conversely, if the system in Fig. 1b was the reference 
with a function H(z), the filter function in Fig. 1a has 
to be H(z™),. 

Equation (4) shows that downsampling a signal 
by a factor of M is achieved by simply keeping ev- 
ery Mth sample and dropping those in between; (5) 
shows that upsampling by a factor of M is achieved 
by inserting M — 1 zeroes between the original sam- 
ples. 


2.3 Polyphase Representation of Systems 
and Signals 


This section describes the polyphase representations 
of signals, summarizing [10] and [2]. 


| oe ee eee, ie) 
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Figure 2: A full period of a digitized cosinus f(n] split into two 
polyphases AP ln 4 [n] and f, Pin n]. By adding up the polyphases the 
original signal is reconstructed. 


A polyphase representation is a parallel represen- 
tation that can be applied to all digital signals. The 
conventional serial form of a system, with one input 
and one output, is referred to as serial representa- 
tion in this document. A generic signal f[n] can be 


split into an arbitrary number M of polyphases pe 


M-1 ( ) 
a SG nl (6) 
A=0 


The exponent (p) indicates that the associated term 
is a polyphase partition of the original signal. 
Any polyphase is characterized as follows, with 
M specifying the number of polyphases and A = 


0,1,2,..., M —1 determining the time delay (phase): 
n+A] if (n+A) mod M=0 
0 else 


The corresponding descriptions in the frequency do- 
main are: 


-> z AFP) (z (8) 


FM) = 7 f[nM+alz-™ (9) 
n=—oo 

Figure 2 shows a polyphase representation of the 
signal f[n] for M = 2. Polyphase representation 
separates the original signal f[n] into M polyphase- 
signals fe each with a delay of 1 in respect to its 
predecessor. Each of the polyphases has to be sam- 
pled up by M, as z™ in (8) indicates. 

The sampling process has a phase (meaning: de- 
lay) of A equivalent to the polyphase’s index ranging 
from 0 to M—1. The process of upsampling com- 
bined with the delay results in the zeroes and shift 
necessary for the polyphases to reconstruct the orig- 
inal signal again when summed up (see (6) and (8)). 


x(n] H(z) yin] 
X(2) Y(z) 


(a) A system in serial representation. 
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(b) The system in polyphase representation with M = 2. 


Figure 3: An FIR-System in serial representation and 
polyphase representation. In polyphase representation, the filter 
taps are intertwined with zeroes, as indicated by zM. Signals 
x[n] and y|n] are identical for both versions. The delay block 
z_! implements the necessary phase of A = 1. 


Likewise, every FIR-System can be split into a 
polyphase representation. This can be conceptual- 
ized by the fact that an FIR-System is, by definition, 
fully characterized by its impulse response, which 
is identical to the filter taps [9]. By splitting the im- 
pulse response into a polyphase representation, the 
system itself can be split. Consequently, a filter in 
serial representation as depicted in Fig. 3a can be 
divided into its polyphase representation and calcu- 
lated in parallel, as clarified in Fig. 3b. As a result of 
this, the original filter’s taps are distributed among 
the polyphase filters and are filled with intermediate 
zeroes, as described in (7) and illustrated by Fig. 2. 
The serial representation is re-entered by summing 
up the polyphases as described in (6) and (8). In an 
implementation, this is achieved through a summa- 
tion point, as clarified in Fig. 3b. 

The signals Y,(z) in Fig. 3b are polyphase-filtered 
subsignals: 


Ya(z) = X(z) H® (2M) (10) 


As these signals are not upsampled, their function 
term contains z instead of z™. 


3 Filter Banks and Channelizers 


A channelizer is a system that splits a spectrum into 
a variable number of smaller, consecutive subbands. 
These subbands are then shifted downwards (i.e., de- 
modulated) to form baseband signals. 


3.1 Fundamentals 


At its core a channelizer can be imagined and imple- 
mented as a cascade of filters cutting out the desired 
spectrum for every output channel. A cascade of 
filters is referred to as a filter bank. Such a system 
is depicted in Fig. 4. A filter bank which slices a 
spectrum into smaller subbands is called an analyti- 
cal filter bank. For simplicity’s sake, filter banks are 
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Figure 4: An analytical filter bank consisting of four channels. 


depicted with only four channels within this docu- 
ment. However, the following derivation does work 
accordingly for larger banks. 


Uniform Filter Banks A uniform filter bank is a fil- 
ter bank with equidistant channels in which all chan- 
nels have the same bandwidth (such as depicted in 
Fig. 5). The sum of all bandwidths amounts to the 
total available bandwidth B = Os: 


M-1 
Bk =B=Os =2n (11) 
k=0 
B 
B. = — 12 
k= (12) 


The channels’ center frequencies Q, are conse- 
quently placed in the middle of their bandwidth 
share: 


fe) 
kM M 


(13) 

To slice the input spectrum continuously, the chan- 
nel borders, and, therefore, the filters’ cutoff frequen- 
cies should ideally be positioned at 


_ ON 


QL. 3 
ae 


(14) 


The transfer functions shown in Fig. 5 adhere to the 
aforementioned rules. 


Filter Banks as a Channelizer’s Groundlayer Fig- 
ure 6 shows a filter bank functioning as a channel- 
izer’s core component. As seen in Fig. 5, the Filters 
H(z) have to be bandpass filters with different cen- 
ter frequencies Ox, retrieving the desired subspectra 
of the input signal X(z). 

Similar to what classical amplitude demodulation 
accomplishes, the demodulation blocks D, are re- 
sponsible for shifting the filtered bandpass signals 
Y, downwards by their center frequency Ox, turn- 
ing them into the baseband signals Y,. In the time 
domain this demodulation is implemented in form 
of a multiplication with a rotating pointer [9]: 


Dy(z) @—0 dg[n] = eft (15) 
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Channel 0 


A(Q) / dB 


-1.00 -0.75 —-0.50 -—0.25 0.00 0.25 0.50 0.75 1.00 
Q/n 


Channel 1 


A(Q) / dB 


-1.00 -0.75 -0.50 —0.25 0.00 0.25 0.50 0.75 1.00 
Q/t 


Channel 2 


A(Q) / dB 


-1.00 -0.75 -0.50 -0.25 0.00 0.25 0.50 0.75 1.00 
Q/n 


Channel 3 


A(Q) / dB 


—1.00 -0.75 —0.50 —0.25 0.00 0.25 0.50 0.75 1.00 
Q/n 


Figure 5: Example of a 4-channel uniform filter bank transfer 
function, implemented by the filter bank of Fig. 4. Each filter 
cuts out a fourth of the available complex bandwidth. 
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Figure 6: The filter bank previously illustrated in Fig. 4, here 
extended to be a channelizer. The demodulation blocks D,(z) 
turn the bandpass signals Y;(z) into baseband signals, which 
then are permissibly sampled down by a factor of four. 


3.2 Downsampling 


The output signals in Fig. 6 each have the bandwidth 
"1, as defined in (12). As these signals are now 
baseband signals, their maximum and minimum fre- 
quencies are 

On 

avi 
Consequently, rate reduction by a factor M is possi- 
ble without violating the Nyquist Theorem: 


(16) 


+Omax = 4 


A=k (7 
While sample rate reductions are widely used in 
practical implementations they are not necessary to 
understand the non-optimized, basic version of a 
Polyphase FFT Channelizer. Therefore, further de- 
tails on rate reductions are provided at a later point 
in the document regarding optimizations in Sec- 
tion 5. 


3.3 Disadvantages of Manual Filter Bank 
Design 


Designing and operating the types of filter banks 
that have been presented so far is difficult because 
of the following reasons: 


¢ M filters have to be designed manually. 

¢ M outputs must be calculated during runtime, 
resulting in the input sequence x[n| having to be 
convolved with all M filter tap sequences M,[n]. 


This paper only considers uniform filter banks. As 
uniform filter banks consist of equally sized chan- 
nels only varying in their center frequencies, it is 
reasonable to dedicate efforts towards reusing a sin- 
gle channel filter. 


3.4 Modulated Filter Banks 


We utilize the properties of the Z-Transform to de- 
rive how shifting a filter’s transfer function H(z) 
corresponds to changes on the filter taps h[n]}: 
H(z)e—oh|n]. By shifting a single Model Filter, all 
other channel filters are obtainable [2], which en- 
ables to design a uniform filter bank. This model 
filter is a lowpass filter that is transformed into band- 
pass filters by shifting it in the frequency domain. 


Modulating Filter Coefficients The modulation 
property! of the Z-Transform is defined as fol- 
lows [12, 9]: 


F (=) e—o a"f[n] VaeCx{0} (18) 


1 Sometimes also referred to as »scaling property« 


If a represents a complex number a = e/® with O 
specifying an arbitrary angle (frequency), the expo- 
nential laws of mathematics show that a becomes a 
complex rotating pointer in the time domain: 


: n 
q” = (ee) = 


One could interpret this utilization of the property 
such that the filter coefficients start »rotating« in the 
complex plane at the frequency ©. By scaling with 
(19), the filter’s transfer function is shifted upwards 
by QO: 


esr (19) 


H (=n) e—oe I hin] = [h[O]e2, h[lJe!2, 
h[2je??, ...] (20) 


Example: Scaling the transfer function Ho(z) of a low- 
pass filter serving as a model filter with a = e7/"/”? 
and operating it at fg = 1000 Hz transforms the low- 
pass into a bandpass with its center frequency at 
250 Hz. 


Uniform Filter Bank Design with Modulated Fil- 
ters The center frequencies of all channels Q, are 
now placed as described by (13). Through (18), a sin- 
gle model filter Hp is transformed into a cascade of 
bandpass filters with complex coefficients. In order 
to design a filter bank with the frequency response 
shown in Fig. 5, the taps for each of the four filters 
must be calculated using those of the model filter as 
a basis. The following equation describes how the 
filter taps for a channel k are computed: 


hyn] = (e#*)" hol (21) 


It is, therefore, possible to create a uniform filter 
bank by modulating a single model filter. 

In the frequency domain representation of (21), we 
define the modulated time discrete complex frequency Z 
in the following form: 


H;(z) = Ho (2 ex ir) =: H(z) (22) 


Designing a uniform filter bank through modu- 
lated filter coefficients has disadvantages, some of 
which already occurred in the conventional filter 
bank: 


¢ It is de facto still necessary to design M filters: 
Each single channel filter’s tap sequence has to 
be computed through (21). 

¢ It is still required to compute M convolutions 
during runtime. 

¢ The originally real-valued filter taps turn into 
complex-valued filter taps, resulting in more 
complicated convolution calculations. 


Nevertheless, modulating a model filter is the first 
step on the way to a more efficient implementation. 


Polyphase FFT Channelizer Derivation 


4 Polyphase FFT Channelizers 


Combining the knowledge from previous sections 
permits designing a highly efficient form of a chan- 
nelizer, referred to as Polyphase FFT Channelizer or, 
briefly, Polyphase Channelizer. 


4.1 Modulation-Independent Filter 
Coefficients 


To derive the new structure, the single channel filter 
H(z) of a channelizer as shown in Fig. 6 is chosen 
with the intention to use it as a universal filter by 
replacing the transfer function by the one described 
in (22). Hence, we utilize a model filter to derive 
all channel filters from it. This filter can now be 
interpreted as a filter responsible for cutting out an 
arbitrary channel, depending on the value of k. With 
adjustments to the transfer function, Fig. 3a could 
represent such a filter. 

The universal filter Hp(Z) is now polyphase- 
decomposed into M polyphases (compare with 
Fig. 3b) according to the rules defined in Section 2.3 
— see especially (8) and (9). As defined previ- 
ously, polyphase decomposition substitutes Z with 
zM_ Therefore, in the modulated universal channel 


2 i ; k M 
filter, Z is expressed as (2 exp(j27rf)) , and ac- 


cording to (8) the universal transfer function in 
polyphase representation for any channel k is 


j= me (2 m, (23) 


The extended form of the summands in (23) is 


; —A M 
z AH) (2M) = (zePuh) “He! ((z | : 
(24) 
The function argument of H®) o, in (24) evaluates to 
2M for every channel k: 


M 
(2 er) = 7M eitk _ 7M for ke No (25) 
Thus, (24) can be simplified to 

eA py) (gM) — (7 o-i2th)\~* py) (Mm 

Zz Hoy (2 Ve (ze Prin) Hoy (2 ) 
=z Hy) cas eb atk (26) 
Equation (26) shows that in polyphase representa- 
tion the coefficients in the modulated filter now are 
independent from any frequency shift and channel in- 


dex k. The coefficients do not rotate anymore as they 
previously did in (20). 


Polyphase FFT Channelizer Derivation 


x[n] on 


(Pp) 4 : (27 OK 
X(2) t >| Ay (2°) >) ed 7 OF 


H(z!) ——» ei 571k 


Y,2) 
j27 2k 


el Mu 


zt 
Zz OE 
$ HY)(24) -——— 
zt a 


L__»] H{P(24) | —__+| of 3h 


Xy 
Figure 7: A modulated universal channel filter in its polyphase 
representation with M = 4. The summation point serves en- 
tering serial representation again. M blocks with X) as inputs 
(highlighted in gray) are necessary to obtain the signals for all 
channels. Except for the value of k, each of the gray substruc- 
tures is identical. 


Consequently, the simplified universal transfer 
function for an individual channel k is: 


NT —ags(P) (pM) (2 Ak 
Halas: )2 Hox (z ) elit (27) 
A=0 


Equations are made more compact to (ultimately) 
obtain a term for the output signal. To do so, the in- 
termediate polyphase-filtered signals X, (see Fig. 7) 
are defined to describe the part of the processed sig- 
nal X which is common across all channels: 


Ma (2\3= Xiz)e-* Ae (2M) (28) 


To obtain the output of a channel k, these intermedi- 
ate signals have to be multiplied with the complex 


value e/i**, Afterwards, the paths are summed up 
to leave the filter’s polyphase representation and re- 
enter serial representation, just like in Fig. 3b. The 
output of an individual channel based on (28) ac- 
cording to (8) then results to: 


M-1 oe 
¥(z) = YY Xa(z) ef a™* (29) 
A=0 


Equations (28) and (29) describe the flowgraph of 
the modulated universal channel filter in Fig. 7. This 
flowgraph illustrates a universal channel filter with 
coefficients that are now independent of k. Chan- 
nel dependency occurs after the polyphase filters 
when the signals X, are multiplied with the channel- 
dependent complex pointers. 

To implement the channelizer depicted in Fig. 6 
using modulated filters, M of the complex pointer 
blocks on the right hand side of X, are necessary. 
These blocks are highlighted in gray in Fig. 7. 

Implementing said pointer blocks manually 
would entail further implementation efforts. Fortu- 
nately, a ready-to-use algorithm implementing the 
desired operations is known since many decades. 


4.2 FFT-Modulated Filter Banks 


As was defined in (29), to retrieve the output signal 
Y,(z) for a channel k, the polyphase-filtered interme- 
diate signals X, have to be multiplied by the com- 


plex number e/ ™* and then be added. Transform- 
ing (29) into the time domain results in the following 
representation: 


M-1 a 
yx[n] = > £,[n] elk with n =0,1,...,00 (30) 
A=0 


Except for the missing factor YW, (30) is identical to 
the computation rule for the Inverse Discrete Fourier 
Transformation (IDFT)? [9]: 


1 Nat Bey 
“= 5 yo XelN (31) 
I=0 


In (31), X; denotes a frequency bin with index I, 
while x; is the corresponding time value after trans- 
forming a discrete spectrum consisting of N bins 
back into a time-domain signal, resulting in N sam- 
ples. The IDFT implements (30): 


¢ X can be interpreted as X, from Figures 7 and 8. 

¢ The sum of all polyphases A in (30) can be inter- 
preted as the sum over all rotation steps / in (31). 

¢ In an instant of time no, a channel output sam- 
ple yx is equal to the output sample x; of the 
IDFT. 


Consequently: 
Ll=A,i=k 
Yk = Xi (32) 
M=N 


The Polyphase FFT Channelizer Considering all 
the previously introduced aspects, the new structure 
in Fig. 8 illustrates an implementation where the 
IDFT block replaces the previously necessary mul- 
tiplications with the complex pointers, as well as 
the summation point to exit polyphase representa- 
tion, both depicted in Fig. 7. The flowgraph includes 
blocks for demodulation and sample rate reduction, 
as discussed in Section 3.1. 

It is important to note that there are differences 
to the earlier presented version of the channelizer 
depicted in Fig. 6. The four filters here are not the 
channel filters. They are the polyphase split of a single 
model filter. As a result of this, the filter functions H 
are indexed with the polyphase index A instead of 
the channel index k from here on. The necessity of 


? The typically used indices n,k have been replaced with i,] to 
prevent confusion with the time index n and channel index k. 


x[n] 
X(z) 


I HP (24) 


‘ad H(z) 


P(e) 


i») Hy(2") 


Polyphase FFT Channelizer Derivation 


yo[m] 
v Y (z¥4) 


i) 


ys[m] 
v Y,(z*") 


ya[m] 
4. |} 


: Y,(2"") 


ya[m] 


v 
Y,(Z*") 


Yi, 


Figure 8: The derived channelizer. The four filters are polyphase partitions of the model filter. The signals Y, are the bandpass filtered 


channels, before they are shifted downwards into baseband. 


this index-change in general is also clarified in Fig. 8 
(X, and Y;). 

After the IDFT, the separate channels can be de- 
modulated and, as their maximum frequency is now 
limited, the sample rate can be reduced. After rate 
reduction, the outputs’ Nyquist Frequencies are all 


In = mn. 


Streamed IDFT Strictly speaking, the filter bank 
as described by (30) operates on mathematical se- 
quences, whereas (Inverse) Discrete Fourier Trans- 
formations operate on data vectors. To obtain the 
output vector Y,,, = [9o,91,---,9m_—i] for a point in 
time no, it is necessary to compute (30) for each k. 
Therefore, this procedure needs to be repeated in in- 
tervals of Ts to obtain an output stream. One could 
interpret the procedure as a »streamed IDFT« where 
the input vector changes with every increment of n: 


%\0| = Es 2 [n] ei wii (33) 
i — N Par l 


The streaming aspect is visualized in Fig. 8. The 
polyphase-branches supply the IDFT with M sam- 
ples. The IDFT’s payload is updated and calcula- 
tions are performed with every increment of n. 


4.3. Derivation Summary 


The following list summarizes the significant steps 
in deriving the polyphase channelizer: 


1. A model filter is modulated with channel index 
k acting as a variable for specifying the target 
channel frequency. (22). 

2. The modulated model filter is split into its 
polyphase representation. (23), (24). 

3. This split results in equal filter coefficients 
across all channels k. (25), (26). 


4. By reshaping the modulated polyphase delays 
z-* from (23) a complex pointer is produced 
that depends on both variables k and A. (26). 

5. This pointer combined with the sum (Fig. 7) nec- 
essary to leave polyphase representation again 
results in equivalence (except for “\) with the 
IDFT-Algorithm. (32). 


It is important to highlight that the IDFT does 
not modulate the filter coefficients — it implements 
the modulation of the polyphase delays 2~* in (23). 
As previously shown through (25), modulating the 
coefficients is no longer necessary, as they became 
independent from any channel index. All of the 
steps listed above, in particular entering and leaving 
polyphase representation, bring forth a mathemati- 
cal expression equivalent to the IDFT. 


4.4 Advantages 


The presented structure provides a variety of ad- 
vantages. For instance, the signals x , are identical 
across all channels and, therefore, only have to be 
calculated once. Furthermore, the IDFT is among 
the most essential mathematical operations in com- 
munication technology. Due to their long history of 
usage, optimized DFT-libraries exist for virtually all 
relevant digital systems. These highly optimized im- 
plementations make use of an algorithm referred to 
as Fast Fourier Transformation (FFT). Available FFT- 
Libraries are likely to be far more efficient than cus- 
tom implementations of the sum described in (30) 
would be. This is why polyphase channelizers are 
predominantly implemented with an (I)FFT. Conse- 
quently, they are often also referred to as Polyphase 
FFT Channelizers. 

Additionally, design effort is reduced as the en- 
gineer has to design only one low pass filter (the 
model filter) and specify the number of channels. At 
this point, however, the total number of filter taps 
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and, consequently, filter calculation effort is not re- 
duced: The model filter’s taps are distributed to the 
polyphases, but have to be filled with zeroes (see 
Fig. 2 and (5)), resulting in the same total number 
of taps. Section 5 shows how these zero-taps can be 
eliminated. 


4.5 IDFT or DFT 


Relevant sources (e.g., Fliege [2] and Harris [6, 8]) al- 
ways present the channelizers with the IDFT. How- 
ever, DFT and IDFT are entirely interchangeable. 
The key to understanding this lies in Section 3.4: It 
was implicitly decided to shift the model filter up- 
wards by setting variable a in (18) to the negative 
pointer e~/k. If the model filter were to be shifted 
downwards with a positive pointer e/k, and the rest 
of the derivation is carried out in the same manner 
as described in previous sections, the result in Fig. 8 
would contain the DFT instead of the IDFT. Due to 
the different shifting direction, this would cause the 
spectral shares of channels one and three in Fig. 5 to 
switch positions. 

The DFT would implement the desired modula- 
tion algorithm exactly, because it lacks the damping 
by a factor “\ which occurs in the IDFT [9]. As this 
difference is not very significant, this paper uses the 
IDFT to maintain congruency with previous publica- 
tions. 

Engineers presumably often prefer the IDFT be- 
cause it suits the following perception better: a chan- 
nelizer slices a bandwidth into channels and de- 
modulates them, so they become baseband signals. 
Hence, one might interpret the IDFT’s input data as 
signals separated in the frequency domain and the 
output as time domain data. Thus, the IDFT acts as 
a separation mechanism between frequency domain 
on the left and the time domain on the right [1] (com- 
pare with Fig. 12). Using a DFT instead would not 
allow for this interpretation to be valid, although the 
system’s behavior would be almost identical. How- 
ever, this interpretation is purely aesthetical and nei- 
ther necessary to design nor to understand the chan- 
nelizer. Section 5.2 provides further details on the 
IDFT’s interpretation. 


5 Polyphase Channelizer 
Optimizations 


This last section introduces the most basic form of an 
optimized channelizer, the »Commutated Polyphase 
Channelizer« and a special version with a double 
commutator that has reduced aliasing compared to 
the primer. 
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(a) Demodulation followed by (b) Rate reduction moved to 
rate reduction. the D-block’s left. 


Figure 9: The rate reduction rule from Section 2.2 applied to the 
demodulation blocks. The D-Block’s transfer function becomes 
equal to 1 and, therefore, obsolete. 
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Figure 10: An intermediate step to the optimized version: the 
rate reduction blocks have first been moved over the D-Blocks, 
making the later obsolete. In a second step, they have been 
moved over the polyphase filters, causing H? Y (z*) to change 
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5.1 Commutated Polyphase Channelizers 


As mentioned in Chap. 1, practical implementa- 
tions typically move the rate reduction to the sys- 
tem’s start to reduce calculation effort. Specialized 
channelizers implementing rate reductions of factors 
M, 5 M and 3 M have been described [13]. As a first 
introduction to the general concept, this section de- 
scribes how the presented channelizer with a rate 
reduction of factor M can be optimized. This sec- 
tion starts off by providing an introduction to the 
general concept of commutated polyphase channel- 
izers and describes the necessary optimization steps 
for the presented channelizer with a rate reduction 
of factor M. A justification is given for why rate re- 
ductions unequal to M are advantageous in the first 
place — afterwards, one such implementation, with a 
rate reduction of a factor ™%, is introduced. 

The following list enumerates the steps necessary 
to optimize the channelizer from Fig. 8. 


1. When the D-Blocks are implemented as rotating 
pointers as in (15), shifting the rate reduction to 
their left results in the pointers becoming equal 
to 1, as only every Mth sample is kept [14]: 


The disappearance of the D-Blocks is visualized 
by Fig. 9. 
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Figure 11: An optimized version of the channelizer. The com- 


mutator moves downwards with a frequency of fs, starting at 
the top again after reaching the last position. 


2. As the IDFT is only a sum and multiplication 
with (per channel k) constant complex numbers, 
the rate reduction can be moved over it without 
further efforts. 

3. When the rate reduction is further moved to 
the polyphase-filters’ left, those filters’ expan- 
sion by factor M disappears, as detailed in Sec- 
tion 2.2: 

H®) (2M) 2, H(z) (35) 
This simplifies computing the filter convolution, 
as coefficients equal to zero disappear and the 
number of taps is reduced by a factor of M (see 
(5) and (4)). 

4. Now, the rate reductions are directly at the sys- 
tem’s start, on the right side of the z~! delays, 
as portrayed in Fig. 10. The rate reductions can 
now be combined with the delays to form a com- 
mutator. 


Ultimately, the steps above result in the Commu- 
tated Polyphase Channelizer depicted in Fig. 11. The 
great advantage of this structure compared to its 
predecessors is a substantial reduction in calculation 
costs: the number of taps is reduced by M. 


Performance Fliege [2] calculates that a channelizer 
such as the one illustrated in Fig. 11, with 32 chan- 
nels and a model filter with 256 taps reduces calcu- 
lation effort by a factor 1046 compared to a structure 
with 32 independent bandpass filters with the same 
number of taps. In general, the efficiency can be at- 
tributed to the following points: 


e Instead of M, only one filter has to be calcu- 
lated during runtime. This filter is split into its 
polyphases. 

¢ The polyphase split allows parallel filter calcula- 
tion, which utilizes resources such as processor 
cores more efficiently. 
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Figure 12: The IDFT can be interpreted to act as a mechanism 
splitting the system into a frequency and a time domain. In their 
respective domains, input and output can be conceptualized as 
vectors. Thus, X and Y are corresponding with each other. 


¢ Moving the rate reduction to the system’s start 
causes 


— the system to be operated at a lower fre- 
quency. 

-— the expansion z™ to collapse to z. Hence, 
there are no more filter taps with zero 
value anymore, which, unfortunately, had 
to be considered when calculating the filter 
convolution. 


¢ The IDFT can be implemented as an IFFT. FFT- 
implementations are available virtually every- 
where and can be trusted to be very performant. 


There seems to be a common consensus that the 
most efficient mechanism to split a signal into an 
even number M of equally large frequency slices is 
a Commutated Polyphase FFT Channelizer. [2, 15, 8, 
16] 


5.2 IDFT as a Gateway to DFT-Properties 


For example when using reduction rates other than 
M, optimized versions may utilize another mathe- 
matical law to move rate reduction-related opera- 
tions to the left side of the IDFT. As visualized in 
Fig. 12, the IDFT permits viewing the channelizer as 
a system separated into frequency domain (on the 
left) and time domain (on the right). Thus, one can 
interpret the transform’s input bins as an input vec- 
tor X and the output samples as an output vector Y. 
As the figure establishes, through the IDFT, the two 
vectors are correspondences to each other: 


Ke—oY (36) 


Therefore, all changes to Y can be replaced by 
changes to X and vice versa, if there is an associated 
property for DFTs. The following section shows how 
the frequency shift property helps transforming rate 
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Figure 13: Symbolic depiction of how a channel’s transfer func- 
tion distorts itself after reducing the sample rate. The gray area 
is the mirrored part of the transfer function exceeding the new 


Nyquist Frequency ft, = 77 fN. 
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(a) Demodulation followed by (b) Rate reduction moved to 
rate reduction. the D-block’s left. 


Figure 14: Moving the rate reduction of “2 as previously in 
Fig. 9 causes the reduction to become a multiplication by +1. 


reduction-related multiplications on Y into bin-swap 
operations on X. This way, an efficient channelizer 
that implements a rate reduction of a factor 4 gets 
constructed [8]. 


5.3 Twice Commutated Channelizer 


The so far presented channelizer has one major dis- 
advantage, which is depicted in Fig. 13: It is critically 
sampled. As filters with infinitely steep flanks are 
impossible to implement, a part of the filter’s trans- 
fer function exceeds the new Nyquist Frequency fy. 
Thus, channels in a channelizer with a rate reduction 
of a factor of M distort themselves. Effectively, this 
leads to aliasing frequencies, being stronger near fy. 


Depending on the available flank steepness, alias- 
ing might or might not be a problem. In any case, if 
increasing the flank steepness is not an option, the 
problem can be lessened through a rate reduction by 
a factor smaller than M. But, as a strong rate reduc- 
tion is desirable due to reduced energy consumption 
and computation cost, it is also not a desirable way 
not to decrease the rate at all. Therefore, a compro- 
mise is adopted: The rate is reduced, but instead of 
choosing the maximum reduction by factor M, the 
second strongest possible reduction by factor ™% is 
applied. The following sections show how such a 
reduction can be implemented efficiently. 


Demodulation Once again, the goal is to eliminate 
the demodulation blocks from Fig. 8 by moving the 
rate reduction over them. Moving the reduction by 
M = ™ over the blocks delivers the following com- 
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Figure 15: Moving the rate reduction to the left turns the D- 
Blocks into multiplications with +1. 
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Figure 16: Vector-based representation of Fig. 15 (compare also 
with Fig. 12). X (following the rate reduction) contains the pay- 
load getting transformed by the IDFT, Y is the output vector. To 
implement the demodulation, elements in Y with uneven chan- 
nel index k must be multiplied with +1. The structure's inputs 
£, are the polyphase-filtered polyphase signals (compare Fig. 8). 


putation rule to get demodulated output samples: 


;27 2 < 
yx[n] = peinje™ => y,[n] = H[Mn] ei" (37) 


with eit = +1. if kn even (38) 
—1 else 


Note that the sign of 7 depends on both k and n. 
Fig. 14 shows how shifting the rate reduction over 
the D-Blocks works in principle. Figure 15 shows 
how shifting rate reductions by a factor ™ to the left 
in the channelizer in Fig. 8 changes the flowgraph. 

Figure 16 and Table 1 clarify Equations (37) and 
(38). As illustrated, the demodulation can be imple- 
mented as a multiplication with —1 for all uneven 
channel indices k and uneven running indices n. 

For the next step, it is helpful to view the IDFT 
as a separator between frequency and time domain 
with input and output vectors, as introduced in Sec- 
tion 5.2. While replacing the D-Blocks with com- 
paratively simple multiplications with +1 is already 
an improvement, there is room for further optimiza- 
tion. The separation in frequency and time domain 
allows utilization of the frequency shift property for 


n 
k 0 1 2 3 
0 1 121 éid31 
11-1 1 =-1 
2 1 1 1 «21 
3 1 -1 1 -!1 


Table 1: Time visualization of Equation (38). Samples of un- 
even channels k are multiplied with —1 for uneven values of n. 
Depicted are the states that the multipliers in Fig. 16 take over 
time. 
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Figure 17: Moving the multiplications in the time domain with 
+1 also to the IFFT’s left is equivalent to switching X’s samples 
(green arrows) in the frequency domain. This structure’s inputs 
and outputs are identical to the ones in Fig. 16. The multipli- 
cation with +1 is equivalent to swapping the upper and lower 
half of X. 


time-discrete Fourier Transformations [9] 
F(wel) e—o e!® Fn]. (39) 


Looking at Figures 16 and 17, one can interpret vec- 
tors X and Y to be corresponding with each other 
through this property. Hereby, X is seen as a vector 
containing frequency bins, and Y as one containing 
samples in time domain. By setting Q to 27 or 7, 
respectively, the complex numbers in the property 
are set to implement the desired multiplication with 
sal 


1= el27 = elOs 


—-l1= el = el 3 Os 


The property permits shifting the multiplication 
from the right side of the IDFT to the left side. 
It shows that multiplying every second sample in 
the output vector with —1 (in the time domain) is 
equal to applying the maximum possible frequency 
change in the input vector. Hence, the frequency 
bins have to be shifted in their vector as far as possi- 
ble. 


Polyphase FFT Channelizer Derivation 


¢——> zen —_»} [2 | —__» 


x[n] 
________-@ 

X(z) 
o—___ >} zm) +—_——| | M/2 ———»> 


e > Zz (M2 +1) ' M/2 > 


1p} z(t) | {we [> 


Figure 18: The rate reductions on the right hand side of the 
polyphase delays. Here presented in absolute numbers and for a 
general bank of width M. The marked delays will be split in the 
next step. 


In summary, the demodulation of a filter bank 
with the sample rate reduced by ™% can be imple- 
mented by shifting the bins of the input vector in 
front of the IDFT. As Fig. 17 illustrates, this results 
in swapping the upper and lower half of the buffer. 
[8] 

Furthermore, Table 1 establishes that the swap op- 
eration in Fig. 17 is only performed for uneven val- 
ues of n. 

Both implementations are in theory equivalent. 
However, it is reasonable to predict higher efficiency 
with the presented method when implemented, i.e., 
in FPGAs: Already available data need only be redi- 
rected into a different memory cell, instead of requir- 
ing a costly multiplication. As digital numbers are 
often implemented using Two’s Complement Encoding, 
it is typically not possible to implement the multipli- 
cation by +1 by merely toggling an algebraic sign 
bit [17]. Therefore, the presented optimization is de- 
sirable. 


Double Commutator In this next step, the sample 
rate reduction is moved to the left-most position, 
similarly as it was done previously in Equation (34) 
and Fig. 10. Moving the reduction by a factor of 
M4 over the polyphase-filters causes the expansion to 
collapse to 


HOM es 2"), 
Finally, the rate reduction is at the left most position, 


on the right-hand side of the polyphase delays. For 
the last optimization step, a different representation 
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Figure 19: Using equivalences to move the rate reduction over 
a part of the delay. 
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Figure 20: Moving the rate reduction by - over the z~? de- 
lays causes the latter to turn into z~' delays. Now, upper and 
lower halves of the flow graph are identical on the rate reduc- 
tion’s left. 


is selected: Figure 18 presents this intermediate state 
of the structure for a bank of width M. The delays 
are presented in absolute numbers instead of build- 
ing them with consecutive z~! delays. The entire 
structure is conceptually split into two halves, with 


their border between z~~! and z~?. The lower 
half of the delays then gets split by extracting z~ 2: 


Fa OE a a ace cet aa (40) 


This delay extraction is visualized in Fig. 19. Once 
more, we apply the rules from Section 2.2: the rate 
reductions in the lower half of Fig. 18 are moved to 
the left. This causes the delays to become equally de- 
laying by a factor of 1. Figure 20 shows the interme- 
diate result of this step. With this figure, the deriva- 
tion is at a very similar point as it was in Fig. 10: it 
is now possible to replace the delays and rate reduc- 
tions with a commutator. In this case, however, two 
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commutators are necessary. The upper and lower 
halves of the structure differ in the lower half’s val- 
ues being delayed by one time step. 

The final result can be seen in Fig. 21, which 
shows the derived structure for M = 4. As men- 
tioned, the buffer switching portrayed in Fig. 17 only 
occurs for uneven n, as Equation (38) demands. Fig- 
ure 22 depicts the system’s output spectra? after the 
system has been stirred up with a Dirac-Impulse. 
Now, the channel borders do not lie at LON any- 


more, but at +9N, Therefore, the self distortion 
portrayed in Fig. 13 has been lessened, while at the 
same time the sample rate is reduced before the fil- 
ters’ costly computations. The filters’ expansion has 
also been reduced to z*. The channelizer presented 
in Fig. 21 is the final result of the optimizations per- 
formed in this document. 


6 Conclusion 


The presented paper revisited Polyphase FFT Chan- 
nelizers roughly fifty years after their first publica- 
tion. It is the author’s wish that this document 
supports other engineers in understanding these ar- 
guably very complicated structures. 
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